A Bose-Einstein Condensate in a Uniform Light-induced Vector Potential 
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We use a two-photon dressing field to create an effective vector gauge potential for Bose-condensed 
*^Rb atoms in the F — 1 hyperfine ground state. The dressed states in this Raman field are spin 
and momentum superpositions, and we adiabatically load the atoms into the lowest energy dressed 
state. The effective Hamiltonian of these neutral atoms is like that of charged particles in a uniform 
magnetic vector potential, whose magnitude is set by the strength and detuning of Raman coupling. 
The spin and momentum decomposition of the dressed states reveals the strength of the effective 
vector potential, and our measurements agree quantitatively with a simple single-particle model. 
While the uniform effective vector potential described here corresponds to zero magnetic field, our 
technique can be extended to non-uniform vector potentials, giving non-zero effective magnetic 
fields. 



Ultracold atoms are an appealing system for the 
study of many-body correlated states relevant to con- 
densed matter physics. These widely tunable systems 
in nearly disorder-free potentials have already realized 
one-dimensional Tonks- Girardeau gases flj, and the su- 
perfluid to Mott-insulator transition of a Bose-Einstein 
condensate (BEG) in an optical lattice The imple- 
mentation of these simple iconic condensed matter sys- 
tems paves the way for more interesting systems with 
exotic correlations and excitations, as in the fractional 
quantum Hall effect (FQHE) of a two-dimensional elec- 
tron system in a strong magnetic field @. Simulating 
such a system with neutral atoms requires an effective 
Lorentz force, which is associated with a vector gauge 
potential. Current experimental approaches involve ro- 
tation of trapped BEGs [1, where low field effects, 
such as the formation of an Abrikosov vortex lattice, have 
been observed. For technical reasons, this approach is 
limited to modest effective fields, too small for FQHE 
physics d, Q. Most recent proposals to create signif- 
icantly larg er effective magnetic fields without rotation 
0, H, 11, [13, Ell involve optical coupling between internal 
atomic states, where the atoms are dressed in a spatially 
dependent manner. The effective magnetic field in such 
light-induced gauge potentials can be understood as a 
consequence of changing into a spatially varying basis of 
internal states. This is similar to the effective magnetic 
field in a rotating BEG arising from changing into the 
frame rotating with the BEG. 

Here we report a first step toward using light-induced 
vector potentials to simulate magnetic fields. By dress- 
ing a BEG with two counter-propagating Raman laser 
beams, we realize a spatially uniform vector gauge po- 
tential. The Raman beams couple internal (spin) states 
with linear momenta differing by twice the photon mo- 
mentum. This gives rise to a spatial gradient (along the 
Raman beam direction) of the phase difference between 
spin components of the dressed state. As we will show, 
this spatially varying state leads to a non-zero vector 
potential when the coupling is detuned from Raman res- 
onance. We adiabatically load the BEG into the dressed 
state, and measure properties of the dressed-state disper- 



sion relation by probing its spin and momentum decom- 
position. Although the atoms are stationary in the lab 
frame (having zero group velocity) , the momenta of the 
individual spin components composing the dressed state 
show a non-zero phase velocity which depends on the 
strength and detuning of the Raman coupling. Our mea- 
surements agree with a simple single-particle model, and 
demonstrate the presence of an effective vector potential. 

We dress a ®^Rb BEG in the F=\ ground state with 
two Raman laser beams counter-propagating along x. 
Together these beams couple states \mp,kx) differing 
in internal angular momentum by h (Am^? = ±1), 
and differing in linear momentum kx by 2kr. Here, 
kr = 27r/A is the single-photon recoil momentum, and 
A is the wavelength of the Raman beams. We de- 
fine Er = h'^k'^/2m as the recoil energy. The family 
of three states coupled by the Raman field, ^'(fca;) — 
|| - l,fc^ -t-2fc^), \0,kx), \ + 2A:,.)|, are labeled by 

the quasi- momentum, fc^;. The Raman beams have fre- 
quencies of lul and (lul + Awl) {Alul > 0), and a bias 
field Boy produces a Zeeman shift fkuz ~ gjisBo ~ 
HAujl (see Fig. [T^b). Since the momentum transfer is 
only along x, the single-particle Hamiltonian can be writ- 
ten as H = Hi(fc^)-h [fi'^ik'^ + kl)/2m + V{r)\ (g)!, where 
Til is the Hamiltonian for the Raman coupling, the Zee- 
man energies and the motion along a;, and 1 is the 3x3 
unit matrix acting on the spin space. V{f) is the state- 
independent trapping potential (arising from a far-off- 
resonance dipole trap and the scalar light shift of the 
Raman beams), and m is the atomic mass. In the rotat- 
ing wave approximation for the frame rotating at Awi, 
Tii/h expressed in the state basis of the family '^{kx) is 
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Here S = {Auj^ — loz) is the detuning from Raman res- 
onance , is the resonant Raman Rabi frequency, and 
e accounts for a small quadratic Zeeman shift (Fig. [TJd) . 
For each k^, diagonalizing TLi gives three energy eigen- 
values Ej(kx) [j = 1, 2, 3). For dressed atoms in state j, 
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(b) Level diagram 
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FIG. 1: (a) The **^Rb BEC in a dipole trap created by two 1550 nm crossed beams in a bias field Boy (gravity is along — z). 
The two Raman laser beams are counter-propagating along x, with frequencies ujl and (ujl + Aoji,), linearly polarized along z 
and y, respectively, (b) Level diagram of Raman coupling within the F = 1 ground state. The linear and quadratic Zeeman 
shifts are uJz and e, and 5 is the Raman detuning, (c) As a function of Raman pulse time, we show the fraction of atoms 
in \mF = — = 0) (open circles), |0, — 2fcr) (solid squares), and | + 1,—Akr) (solid triangles), the states comprising the 
'^{kx ~ —2kr) family. The atoms start in | — l,fca; = 0), and are nearly resonant for the | — 1,0) ~+ |0, — 2fcr) transition at 
hS — ~4.22Er. We determine hQa = 6.63(4)i5r- by a global fit (solid lines) to the populations in ^(— 2fcr). 



Ej{kx) is the effective dispersion relation, which depends 
on experimental parameters, (5, VLr, and e (left panels 
of Fig. (2). For example, the number of energy minima 
(from one to three) and their positions, fcmin, are experi- 
mentally tunable. Around each ftmin, the dispersion can 
be expanded as E{kx) ~ f^{kx — fcmin)^/2m*, where m* 
is an effective mass. In this expansion, we identify fcmin 
with the light-induced vector gauge potential, in analogy 
to the Hamiltonian for a particle of charge q in the usual 
magnetic vector potential A: [p — qA)'^/2m. In our ex- 
periment, we load a trapped BEC into the lowest energy 
j = 1 dressed state, and measure its quasi- momentum, 
equal to fcmin for adiabatic loading. 

Our experiment starts with a 3D *^Rb BEC in a 
combined magnetic-quadrupole plus optical trap, as de- 
scribed in Ref. [l^. We transfer the atoms to an all- 
optical crossed dipole trap, formed by a pair of 1550 nm 
beams, by ramping the quadrupole field gradient to zero. 
The trapping beams are aligned along x — y (horizontal 
beam) and at ~ 10° from z (vertical beam). A uniform 
bias field along y gives a linear Zeeman shift loz/^t: ~ 
3.25 MHz and a quadratic shift e/27r = 1.55 kHz. The 
optically-trapped condensate typically has N = 2.5 x 10^ 
atoms in \mp — —l^k^ = 0), with measured trap fre- 
quencies of « 30 Hz parallel to, and « 95 Hz perpendic- 
ular to, the horizontal beam. 

In order to Raman-couple states differing in mp by 
±1, the A — 804.3 nm Raman beams are linearly polar- 
ized along y and i, corresponding to tt and cr relative 
to the quantization axis y. The beams have 1/e^ radii 
of 180(20) ^m 15], much larger than the 20 ^im BEC. 
These beams give a total scalar light shift up to 60Er, 
where Er = h x 3.55 kHz, and contribute an additional 
harmonic potential with frequency up to 50 Hz along y 
and z. The differential light shift between adjacent mp 
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FIG. 2: Left panels: Energy-momentum dispersion curves 
E(kx) for he = 0.44£'r and detuning ft5 = in (a) and h5 — 
— 2Er in (b). The thin solid curves denote the three states 
I — 1, kx + 2kr), |0, fca;), I + 1, kx — 2kr) in absence of Raman cou- 
pling; the thick solid, dotted and dash-dotted curves indicate 
dressed states at Raman Rabi frequency Qr = 4.85Er/h. The 
arrows indicate kx = fcmin in the j = 1 dressed state. Right 
panels: Time-of-fiight images of the Raman-dressed state at 
hQ.R = 4.85(35)£;r., for detuning M = in (a), and hS = -2Er 
in (b). The Raman beams are along x, and the three spin 
and momentum components, | — 1, fcmin + 2fcr), |0, fcmin) and 
I + 1, fcmin — 2fcr), are separated along y (after a small shear 
in the image which realigns the Stern- Gerlarch gradient di- 
rection along y). 
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FIG. 3: Quasi-momentum kx of the Raman-dressed state ver- 
sus sweeping time taw of Raman detuning from fij = to 
hS = (solid circles) and hS = 2Er (open squares), at 

hQa = 5Er. The crosses denote kx = without sweeping 
{hS = 0). The inset shows the time sequence of and S for 
the loading into a detuned Raman-dressed state. 



states arising from the combination of misalignment and 
imperfect polarization is estimated to be smaller than 

0. 2Er. We determine fi/j by observing population oscil- 
lations driven by the Raman beams and fitting to the 
expected behavior, as shown in Fig. [It. 

We developed a procedure to adiabatically load the 
\ — l,kj. = 0) BEC into the lowest energy, j = 1, Raman- 
dressed state. For a BEC initially in |0,fcj, = 0), this 
could be achieved simply by slowly turning on the Ra- 
man beams at detuning S — 0, resulting in the j — 

1, kx = dressed state, located at the minimum of 
Ej=i{kx). However, our initial state | — l,fc2, = 0), for 
which kx = kx + 2kr, is a member of the family of states 
"^{kx — —2kr). This simple loading procedure would 
therefore result in the j = 1, kx — —2kr dressed state. 
In order to transfer our kx — —2kr BEC to the desired 
kx = dressed state, we use an additional rf coupling 
which mixes families of states with kx differing by 2kr- 
The rf frequency is always equal to the frequency differ- 
ence between two Raman beams, Aujl/27t — 3.250 MHz. 
The loading sequence is as follows: (i) We turn on the rf 
coupling in 1 ms to a resonant Rabi frequency Qj:{/2tt— 
12 kHz, with an initial detuning H6 = 15Er. We then 
sweep the detuning to resonance in 9 ms {S is var- 
ied by ramping the bias field Bq, leaving Alol constant.), 
loading the atoms into the lowest energy rf-dressed state, 
a spin superposition still at fc^; = 0. (ii) We ramp on the 
Raman coupling in ton = 20 ms to a variable Rabi fre- 
quency fifl, and then turn off the rf in 2 ms. This loads 
the BEC into the j = 1, fc^, = dressed state. In the 
dressed state, the BEC heats at ^ 0.3 nK/ms, we believe 
due to technical noise in Bq or Alol. In addition, we 
observe unwanted population in the j = 2 dressed state. 



Therefore, as part of step (ii), we further evaporatively 
cool by decreasing the intensity of the horizontal trapping 
beam by 20% in 20 ms, after which it remains constant. 
The evaporation typically reduces the BEC number to 
7 X 10". 

In the final step (iii), we transfer the fc^: = dressed 
state into a non-zero kx state by sweeping the detuning 
from resonance to a variable S in tsw and holding for 
th. As we show below, for tg^ > 20 ms, this process 
adiabatically loads the atoms into the j = 1 Raman- 
dressed state at kx = kminiS)- 

We characterize the Raman-dressed state by abruptly 
turning off the dipole trap and the Raman beams in less 
than 1 fis, projecting the atomic state onto its individ- 
ual spin and momentum components. The atoms then 
expand in a magnetic field gradient applied during time- 
of-flight (TOF) approximately along y, and the three spin 
states spatially separate due to the Stern-Gerlach effect. 
Imaging the atoms after a 20 ms TOF gives the momen- 
tum and spin composition of the dressed state (see right 
panels of Fig. [J). The quasi-momentum kx of the BEC 
is given by the shift along x of all three spin compo- 
nents from their positions when fc^ = 0. We determine 
kx = to within 0.03 kr by using the loading procedure 
above, except that in (ii), the Raman frequencies are suf- 
ficiently far-detuned that the atoms experience the scalar 
light shift but are not Raman-coupled; then the rf and 
Raman fields are snapped off concurrently before TOF. 

It is straightforward to adiabatically load into the rf- 
dressed state in step (i). To ensure the subsequent adia- 
batic loading from the rf-dressed state into the Raman- 
drcsscd state in step (ii), we ramped up the Raman power 
in a time ton, then ramped it off in ton- For ton > 2 ms, 
there is no discernable excitation into other rf-dressed 
states. We conservatively use ion = 20 ms as the loading 
time in step (ii). To determine tsw in step (iii), we mea- 
sure kx as a function of tgw for negligible t^. Figure [3] 
shows kx versus tg-w at M — 0, 1, 2Er and hflu = 5Er, 
for th — 0.1 ms. We use tsw — 20 ms, the time at which 
kx has nearly reached its equilibrium value. This adia- 
batic following is enabled by the external trap, and the 
time scale is comparable to our typical trapping periods. 
After tsw, we add a hold time th = 20 ms during which 
residual excitations damp. 

Figure [2] shows spin-resolved TOF images of adiabati- 
cally loaded Raman-dressed states at hilj^ ~ 4.85(35)£'r 
for hS = —2 and OEr- The resonance condition, (5 = 0, 
is determined from the symmetry of the rf-dressed state 
with an uncertainty of h x 1.5 kHz=0.4£'r, limited 
by the stability of Bq. The quasi-momentum kx = k^in, 
measured as a function of S, is shown in Fig. |4^, along 
with the calculated fc,iiin(^)- Each different fcmin repre- 
sents a spatially uniform vector gauge potential, analo- 
gous to the magnetic vector potential, with /cmin = qA/h. 
For such a uniform A, the magnetic field B = W x A 
would be zero. 

For Rabi frequency flji < AATEr/h, Ej=i{kx) has mul- 
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FIG. 4: (a) Measured quasi-momentum kmin versus detun- 
ing 6, and tlie calculated fcmin at Rabi frequency Qr — 
4.85(35)iJr/ft (with the uncertainty indicated by the shaded 
region). The inset shows the time sequence of and S. 
(b) Quasi-momentum fcmini (circles) and fcmin2 (triangles) 
versus detuning 5, and the calculated fcmin (solid lines) at 
hilu = 2.95(35)i5r. The intensity of the circles indicates the 
population of |0, fcmini) relative to its value ai S — 0. The 
intensity of the triangles for S < {5 > 0) indicates the pop- 
ulation of I -I- l,fcmin2 — 2kr} (| — l,fcinin2 + 2fc,.)) relative to 
its value at hS = —3Er {hS = 3Er). The insets picture E{kx) 
at hS = —0.8 and —2.5Er as indicated by the arrows. 



tiple minima for some detunings 6. For data taken as 
in Fig. [Da,, Fig. [3J3 shows such multiple fcmin, expected 
for Sa < \S\ < Sb, where HSa = 0.35(20)-Br- and Mb — 
1.2{l)Er, at = 2.95{35)Er. At small and large de- 
tuning, we observe only one k^in- However, at intermedi- 
ate detuning, we observe two fcmin at fcmini and fcmin2 p^ - 
In the figure, the relative population is indicated by the 
gray scale, and the solid line is the calculated value(s). 
The observed values of fcmin(^) agree well with the theory, 
and a simultaneous fit of fcmini and fcmin2 with Qfi as a 
free parameter gives hilfi = 2.9(3)i?r, which agrees well 
with the independently measured hTlfj — 2.95(35)£'r- 
As expected, we observe population in only one for 
(5 1 < (5a, however, for some \d\ > 6b^ there is discernable 
population in two k^ where we expect only one fcmin. We 
attribute this to non-adiabicity when the minimum at 
fcmini disappears during the detuning sweep in loading 
step (iii). 

In conclusion, we have prepared a Bose-Einstcin con- 
densate in a Raman-dressed state with a non-zero quasi- 
momentum, controlled by the Rabi frequency and the 
detuning 5. This technique, in conjunction with a spatial 
gradient of the bias field (therefore the detuning) along 
gives a spatial gradient of light-induced momentum 
fcmin(j/)£ and creates an effective magnetic field along z 
The analog of the magnetic length Ib — \/h/{qB) 
(the classical cyclotron radius of an orbit with one unit 
of angular momentum H) is (9fcmin/9y)~^/^, and is about 
the spacing between vortices in a vortex lattice like that 
formed in a slowly rotating BEC [5]. For a gradient in 
fcmin of one kr across a condensate of radius R, this gives 
an analog magnetic length « 1.4 fim for a i? « 10 /im 
3D BEC, and we expect a ~ 25-vortex lattice to form in 
the condensate. For dilute 2D systems, this is sufficient 
to reach the quantum Hall regime. 
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